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SUMMARY 

A simple model is given for the frequency response of boiler inlet pressure to inlet 
flow. The boiler analyzed is a forced flow, single tube, shell and tube boiler with no 
inserts. The model yields a simple description of the boiler dynamic characteristics. 

The characteristics included in the model are boiler pressure drop, time required for the 
liquid temperature to reach saturation temperature, inertia of the liquid in the subcooled 
region, and mass storage due to the movement of the boundary between the subcooled 
and two-phase regions. The model is compared with experimental boiler frequency 
response data. The boiling fluid was Freon-113, which was heated by counterflowing 
hot water. The boiler exit quality ranged from 12 to 63 percent. Good agreement be- 
tween the model and the data was obtained for qualities ranging from 23 to 63 percent. 


INTRODUCTION 

Two -phase flow stability problems have been encountered in the development of 
boiler systems in which liquid is boiled in heated passages. Boiler instabilities can 
produce degraded system performance and in some cases can lead to system failure. 
Boiler instabilities, which have been detected by many of the investigators listed in 
references 1 and 2, are complex in nature and difficult to analyze. An important class 
of instabilities (manifested by unsteady flow) are those which result from a pressure-flow 
coupling of the boiler and other loop components (ref. 3). To describe this coupling, 
the dynamic characteristics of all the loop components must be known. Dynamic 
techniques have been developed to describe the all-liquid and all-vapor portions of the 
loop (ref. 4). But due to the complexity of two-phase flow (see ref. 5) and the large 
number of possible boiler configurations, a general description of boiler dynamic char- 
acteristics is not available. 



Insight into the dynamic characteristics of boilers can be obtained from the boiler 
frequency response, that is, the response of one boiler variable to small amplitude 
sinusoidal perturbations in another variable. This report presents a simple model for 
the frequency response of boiler inlet pressure to boiler inlet flow. The model includes 
the effects of boiler pressure drop, the time required for the liquid temperature to 
reach saturation temperature, inertia of the liquid in the subcooled region, and mass 
storage due to the movement of the interface between the liquid and two-phase regions of 
the boiler. The model is compared with the experimental frequency response data 
reported in reference 6. 

The data of reference 6 were obtained by imposing small amplitude sinusoidal 
perturbations on the pressure and flow at the inlet of a hollow, single tube, shell and 
tube boiler. A schematic of this boiler is shown in figure 1. Boiler dimensions are 
given in this figure. The boiling fluid, Freon-113, flowed vertically upward and was 
heated by counterflowing hot water. The Freon exhausted into a constant pressure 
plenum tank. Frequency response data were obtained over a frequency range from 0. 04 
to 4. 0 cps (Hz). The data were taken at three nominal values of inlet water temperature: 
210°, 230°, and 250° F (372°, 383°, and 394° K). For each inlet water temperature, 
data were obtained at several Freon flow rates ranging from 175 to 656 pounds per hour 
(0. 022 to 0. 083 kG/sec). The Freon inlet temperature was about 70° F (285° K) for all 
runs. The Freon vapor exit quality ranged from 12 to 63 percent. A detailed descrip- 
tion of the experiment is presented in reference 6. 

ANALYSIS 

The analysis of the boiler dynamic characteristics is based on the following major 
assumptions: 

(1) All of the steady-state boiler pressure drop can be lumped at the boiler outlet. 

(2) The axial density distribution can be approximated by assuming constant densities 
in both the all-liquid and two-phase regions with a discontinuity at the boundary between 
these regions. 

(3) The density in the two-phase region is much smaller than the liquid density. 

(4) The boundary between the all -liquid and two -phase regions is the point where the 
bulk liquid temperature reaches saturation. 

(5) Changes in the position of the boundary between the all-liquid and two -phase 
regions are due only to changes in the flow at the boiler inlet and changes in pressure at 
the liquid two -phase boundary. 

It was also assumed that variations in boiler variables are small amplitude sine waves. 
Thus equations can be linearized and operator notation can be used for time derivatives. 
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A schematic of the boiler indicating the major assumptions made in the analysis is 
shown in figure 2 . 

Liquid flows into the boiler and is heated until its bulk temperature reaches satura- 
tion temperature at x = L . (All symbols are defined in the appendix. ) The saturation 

SC 

temperature is determined by the pressure at the liquid two -phase boundary P T . At 

^sc 

steady state, is assumed to equal the pressure at the boiler inlet P- n , but due to 

SC 

the inertia of the liquid P L will be different from P in when the liquid is accelerated. 

SC 

When the liquid temperature reaches saturation temperature, it begins to boil and its 
density is assumed to drop from P L to P-pp. If the flow into the boiler or the pressure 
at the liquid two-phase boundary changes, the position at which the liquid temperature 
reaches saturation will move. Neglecting density changes in the two-phase region, the 
total amount of mass in the boiler will change by an amount (p^ - P T p)Aj. AL gc . The 
equation for conservation of mass, with P^,p « p^, is written in small perturbation 
form as 

d AL 

AW out “ AW in ' p lA (1) 

From the assumption that all of the steady-state pressure drop can be lumped at the 
boiler outlet, and since the steady-state value of P T is assumed to equal the steady- 
state value of Pj n , it follows that 


AP 


sc 


d < P in - P out> 
dWi. 


AW out = * *W Qut 


Combining equations (1) and (2) results in 


( 2 ) 


AP, 


/ d AL 0 „ 

= R[ AW-. A SC 


SC 


\ 


in 




dt 


(3) 


From the assumption that changes in the location of the boundary between the liquid 
and two-phase regions are functions only of changes in the flow at the boiler inlet and 
changes in the pressure at the boundary, the following is determined: 


AL 


sc 



(4) 
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1 


is the transfer function relating changes in the subcooled 


is the transfer function which 


where (AL „/AW 4 „ 
\ sc' in 


length to changes in boiler inlet flow and ( AL gc / AP^ 


relates changes in subcooled length to pressure changes at the liquid two -phase boundary. 

The transfer function (AL. /AW. ) was derived from the energy equation for the 

\ sc 1: n /p 

^sc 

subcooled region in reference 7. This result is given in the following equation: 


sc / 1 - e 


where T gc = P L A t L sc /W f 


The derivation in reference 7 assumed that the heat flux in the subcooled region was 
not a function of position or time. The steady-state heating fluid temperature profiles 
presented in reference 6 have a slope which is nearly constant in the subcooled region, 
thus indicating that the heat flux did not vary much with position in the subcooled region. 
It is also assumed herein that the heat flux is not a function of time. Thus the derivation 
in reference 7 can be applied to the data of reference 6 and used in this model. 

If it is assumed that the change in boundary position due to a pressure change is not 
a function of frequency, the transfer function LAL^/AP^ \ can be evaluated from 


steady -state data. Thus (aL^/AP^ j = (0L gc /3P L j . Equation (4) can then be 


written as 


sc / 1 - e 


AW in + 


Putting equation (3) in operator notation and combining it with equation (6) give 


Wl-e 8C \ / 3L sc 

AP, = rJaW- - Pt A.JW — ] AW- + SC 

L SC I in L 1 W t l j«r / in 0P T 
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Since V = P L A t L sc /W t , 


AP, 


Re 


-JWT 


SC 


sc 


1 + jwRP L AJ 


'3L 

3Pi 


( 8 ) 


sc 


sc' 


W + 


In small perturbation form, the momentum equation for the all-liquid region is 


AP. = AP, + iui AW. 
in ^sc J Aj. in 


( 9 ) 


Combining equations (8) and (9) yields 


AP 


in 


Re 


-JCOT 


sc 


AW. 

in 


l + jcoRp L A t 


3L 


+ 


sc 


sc 

\ 


3Pt 


sc- 


W. 


(10) 


Equation (10) gives the dynamic ratio of pressure perturbation at the boiler inlet to flow 
perturbation at the boiler inlet. Since this ratio was measured in the experiment 
reported in reference 6, a direct comparison of the model with experimental data can be 
made. 


EVALUATION OF MODEL PARAMETERS 

Since variations in boiler variables were assumed to be small, all of the parameters 
in equation (10) can be defined in terms of steady -state data. The boiler resistance R 
is defined as the slope of the steady -state boiler pressure drop against flow curve. 

Plots of boiler pressure drop against flow, for the boiler used in reference 6, are shown 
in figure 3. Each curve in figure 3 is for constant inlet water temperature. Thus for a 
given inlet water temperature and Freon flow, the boiler resistance can be determined 
from figure 3. 

The subcooled length L is defined as the distance from the boiler inlet at which 
the Freon temperature reaches saturation temperature. This distance can be determined 


5 



from a heat balance between the Freon and the heating fluid (water) and the heating fluid 
temperature profile. Figure 4 is a plot of subcooled length against flow for various 
constant inlet water temperatures. The value of L_„ at a particular Freon flow and 
water temperature was taken from the curve. This value was then used to evaluate r gc , 

T sc = L sc/ V in’ f nd l ) e lic ^ d inertance L sc /A r 

The term (3L gc fop^ J is determined from the assumption that changes in the 


subcooled length are functions of flow and pressure changes only. Thus 


't /P 


Since P L = P in at steady state and P Qut is constant, 
sc 


L sc d < p i„ - p „ut> 


At steady state, ^9L gc /0W t j = (AL gc /AW t j . Thus from equation (5), 


J sc, u>=0 


dL L 
sc _ sc 

dW t W t 


where dL gc /dW^. is obtained from the curves shown in figure 4. Boiler steady-state 
conditions and model parameters are listed in table I. 
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DISCUSSION OF BOILER MODEL AND COMPARISON WITH DATA 


The ratio of boiler inlet pressure perturbation to boiler inlet flow perturbation has 
been derived herein by solving conservation equations for the boiling process. These 
equations were based on a simplified picture of the boiling process. The final result is 
given in equation (10). This ratio of boiler inlet pressure perturbation to boiler inlet 
flow perturbation is referred to as the boiler inlet impedance. The term R in equa- 
tion (10) is referred to as the boiler resistance and is defined as the slope of the steady- 
state boiler pressure drop against flow curve. Since the steady -state boiler pressure 
drop was placed at the boiler exit, pressure perturbations due to this pressure drop 
are produced by flow perturbations out of the boiler. The quantity 


-JWT C 


1 + 


jujRp^Aj. 


8L sc/ 3P L 


SC 


W. 


t-l 


, which multiplies the boiler resistance, is the 


ratio of flow perturbation at the boiler exit to flow perturbation at the boiler inlet. Thus, 
the first quantity in equation (10) is the ratio of the pressure perturbation caused by exit 


-JU)T 

flow perturbation to the inlet flow perturbation. The term e s is the result of mass 

storage in the boiler due to the dependency of the position of the liquid two-phase 

boundary on boiler inlet flow. This term indicates that there is a dead time t be- 

s c 

tween flow out of the boiler and flow into the boiler. That is, the flow out of the boiler 
at time t equals the flow into the boiler at time t - r . For sinusoidal flow 

SC 

perturbations, this dead time produces a phase lag between flow perturbations out of the 
boiler and flow perturbations into the boiler which increases linearly with frequency. The 
dead time has no effect on the amplitude ratio of the flow perturbations. 

The term Rp^A^. ^ is due to mass storage as a result of the 

sc w t 


dependency of the position of the liquid two -phase boundary on pressure. This mass 
storage produces a decrease in the amplitude ratio of the flow perturbations as frequency 
increases and produces up to SO 0 phase lag between flow perturbation out of the boiler 
and flow perturbation into the boiler. The decrease in amplitude ratio is due to the fact 
that, as frequency increases, more of the inlet flow perturbation goes into mass storage 
and less out of the boiler. 

The quantity jajL gc /Aj. comes from the momentum equation for the liquid region of 
the boiler. This quantity has very little effect on the boiler impedance for frequencies 
below 1 cps (Hz) since its magnitude is much less than that of the quantity 

At higher frequencies, where the magnitudes 


Re 


sc 


1 + ^■Rp h A t [dh sc /dP h 


sc 


W, 


t-l 


of these two quantities are nearly equal, the two quantities can subtract or add to produce 
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sharp dips or peaks in the magnitude of the boiler impedance. 

A direct comparison of the model derived in this report with the data from refer- 
ence 6 is presented in figure 5. This figure consists of polar plots of the complex ratio 
of boiler inlet pressure perturbation to boiler inlet flow perturbation. The solid curve 
was obtained from the model (eq. (10)), and the dashed curve was obtained by drawing 
a smooth curve through the data of reference 6. Figure 5 indicates that both the model 
and the data have very similar characteristics over most of the frequency range for these 
nine runs. Two runs not compared (runs 1 and 2 of ref. 6) had low vapor exit qualities 
(12 and 13 percent). The 13 percent run was on the negative slope region of the boiler 
pressure drop against flow curve. The 12 percent run was near this region. Satisfactory 
agreement between the model and the data was not obtained for these two runs. Thus the 
model does not appear to work near the negative slope region of boiler pressure drop 
against flow. Stenning and Veziroglu (ref. 8) also found that a single model would not 
describe boiler dynamic characteristics in both the positive and negative slope regions of 
the boiler pressure drop against flow curve. 

In addition to the comparison in polar form, the model and the data are compared 
in plots of amplitude and phase against frequency in figure 6 for several of the runs. The 
solid curve was obtained from the model. The circles are the actual data points. And 
the dashed curve is a fit of the lower frequency data. It was found that a curve of the 
-jWTj 

form Ke /(I + jcuTg) fit the lower frequency data quite well for all but two of the 

runs (runs 6 and 9) that are compared with the model. The form of the equation used for 

the data fit is identical to the equation derived in this report (eq. (10)) for cases where 

the term jwL^/A^ in the model is negligible. Thus for low frequency, where jwL gc /A t 

is small, the model and the data fit should coincide and there should be a one-to-one 

correspondence between the parameters of the model and those of the data fit. Thus R 

should equal K, r gc should equal Tj, and RP-r A^. (3L SC /3P^ j should equal T 2 if 

\ sc/^ 

the parameters of the boiler model were evaluated properly. Values of r ^ and were 
obtained for the seven runs which could be fit by the equation presented before. Since the 
value of K is determined only by the very low frequency amplitude, K is compared 
with R for all nine runs. Figures 7 to 9 are plots of the parameters in the data fit 
against the corresponding parameters in the model. Good agreement is obtained for all 
three parameters compared. 


RESUME 


The model presented in this report yields a simple expression relating perturbations 
in boiler inlet pressure to perturbations in boiler inlet flow. The model includes the 
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effects of boiler pressure drop, time required for the liquid temperature to reach satu- 
ration, inertia of the liquid in the subcooled region, and mass storage due to the move- 
ment of the boundary between the subcooled and two -phase regions. The model is 
compared with frequency response data reported in reference 6. The model and the 
data have very similar characteristics over most of the frequency range for nine of the 
eleven runs reported in reference 6. In addition, three of the four model parameters 
are compared with corresponding parameters from a fit of the data of reference 6. Good 
agreement is obtained. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 11, 1967, 

120-27-04-27-22. 
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APPENDIX - SYMBOLS 


Aj. cross-sectional area of boiler 

tube, ft 2 (m 2 ) 

j iTi 


L sc 

subcooled length, ft (m) 

AL sc 

perturbation in subcooled 
length, ft (m) 

P in 

pressure at boiler inlet, 
lb/ft 2 (N/m 2 ) 

AP in 

perturbation in boiler inlet 
pressure, lb/ft 2 (N/m 2 ) 

%c 

pressure at x = L , 
lb/ft 2 (N/m 2 ) 

ap l 

SC 

perturbation in pressure at 
x = L gc , lb/ft 2 (N/m 2 ) 

P out 

pressure at boiler outlet, 
lb/ft 2 (N/m 2 ) 

R 

boiler resistance, 
(lb/ft 2 )/(slug/sec) 
[(N/ m 2 )/ (kg/ sec)] 


t 

time, sec 

v in 

liquid velocity at boiler inlet, 
ft/sec (m/sec) 

AW in 

perturbation in flow at boiler 
inlet, slug/sec (kg/sec) 

AW out 

perturbation in flow at boiler 
outlet, slug/sec (kg/sec) 

W t 

boiler flow rate, 
slug/sec (kg/sec) 

X 

position, ft (m) 

P L 

liquid density, slug/ft 2 (kg/m 2 ) 

P TP 

two -phase density, slug/ft 
(kg/m 3 ) 

r sc 

time required for liquid tem- 
perature to reach saturation, 
sec 

C 0 

frequency, rad/sec 
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TABLE I. - BOILER STEADY-STATE CONDITIONS AND VALUES OF BOILER MODEL PARAMETERS 


(a) English units. 


Run 

number 
of refer- 
ence 6 

Freon 

flow 

rate, 

slug/sec 

Vapor 
quality 
at boiler 
exit, 
percent 

Subcooled 

length, 

L sc’ 

ft 

Boiler 

resistance, 

R = dPin/d w t 
(lb/ ft^)/ (slug/ sec) 

Subcooled 
dead time, 

r sc ~ L sc / V in' 
sec 

Inertance 
of liquid 
in sub- 
cooled 
region, 

L sc/ A f 

ft' 1 

Change in sub- 
cooled length 
due to pres- 
sure change, 

( ?L - /ap O ■ 

sc/ w t 

(ft)/(ib/ft 2 ) 

Mass storage 
time constant, 

\ sc/ Wt 

sec 





Nominal water inlet temperature, 210° 



3 

0.00388 

29 

1.92 

1.25x10 s 

0.788 

3560 

2.61X10" 3 

0.518 

4 

.00250 

50 

.84 

4.25 

.535 

1560 

.81 

.547 

5 

.00154 

63 

.26 

6.25 

.271 

480 

.69 

.685 





Nominal water inlet temperature, 230° 

7 


6 

0.00561 1 

23 

2.48 1 

5 

0.85X10 3 

0.702 

4590 

1.38X10 -3 

0. 187 

7 

. 00 395 

37 

1.34 

3.50 

. 542 

2480 

1.20 

.670 

8 

.00250 

45 

.40 

4.75 

.255 

740 

.85 

.640 

9 

.00153 

47 

.10 

3. 50 

. 105 

185 

.21 

. 118 





Nominal water inlet temperature, 250° 

F 


10 

0.00395 

31 

0.60 

4.25X10 5 

0.242 

1110 

1.25X10' 3 

0.842 

11 

. 00520 

31 

1.78 

5.25 

.546 

3300 

1. 18 

.985 


(b) SI units. 


Run 

number 
of refer- 
ence 6 

Freon 

flow 

rate, 

kg/sec 

Vapor 

quality 

at boiler 
exit, 

percent 

Subcooled 

length, 

L sc’ 

M 

Boiler 

resistance, 

R = dP in /dW t , 
(N/ nr)/ (kg/ sec) 

Subcooled 
dead time, 

T sc = ^sc^in’ 
sec 

Inertance 
of liquid 
in sub- 
cooled 
region, 

L sc /A t’ 

M' 1 

Change in sub- 
cooled length 
due to pres- 
sure change, 

sc w t 

m/(N/m^) 

Mass storage 
time constant, 

sc w t 

sec 

Nominal water inlet temperature, 372° K 

3 

0.0566 

29 

0.58 

4. 10X10 5 

0.788 

11 700 

16.6xl0“ 6 

0.518 

4 

.0364 

50 

.26 

13.9 

.535 

5 120 

5. 16 

.547 

5 

.0224 

63 

.08 

20.5 

.271 

1 570 

4.40 

.685 

Nominal water inlet temperature, 383° K 

6 

0.0818 

23 

0.76 

2.78x10 s 

0.702 

15 000 

8. 80X1 O’ 6 

0. 187 

7 

.0576 

37 

.41 

11. 5 

. 542 

8 140 

7.65 

.670 

8 

.0364 

45 

. 12 

15.6 

.255 

2 430 

5.42 

.640 

9 

.0223 

47 

.04 

11.5 

. 105 

607 

1. 34 

.118 

Nominal water inlet temperature, 394° K 

10 

0.0576 

31 

0. 18 

13.9X10 5 

0.242 

3 640 

7.97X10' 6 

0.842 

11 

.0759 

31 

.54 

17.2 

.546 

10 800 

7.53 

.985 
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Figure 1. - Schematic of boiler used in reference 6. 



Boiler pressure drop, N/m 





Subcooled length, ft 



•° 2 -04 .06 .08 .10 

Flow rate, kg/sec 


Figure 4. - Subcooled length as function of flow rate for constant water 
inlet temperature. 
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(b) Run 4. 


Figure 5. - Comparison of boiler model with experimental data. (Steady-state conditions and 
values of model parameters listed in table I. ) 






































Amplitude ratio, | AP in /AW in | . (N/'m) /(kg/sec) 


Q Data of ref. 6 

From boiler model 

— Fit of lower frequency data 



Frequency, eps (Hz) 


Figure 6. - Comparison of experimental data, boiler model, and data fit. (Steady-state conditions and values of model 
parameters listed in table I. ) 
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Figure 7. - Comparison of boiler resistance R (from boiler model) with K (from 
data fit). 



Tg (from data fit), sec 


Figure 9. 
with t. 


- Comparison of Rp L A t (?L sc /rp L ) from boiler model 
(from data fit). ' sc /W t 
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